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Abstract 

We consider a magnetic impurity which interacts by hybridization with a 

system of weakly correlated electrons and determine the energy of the ground 

state by means of an 1/iV/ expansion. The correlations among the conduction 

electrons are described by a Hubbard Hamiltonian and are treated to lowest 

order in the interaction strength. We find that their effect on the Kondo 

temperature, Tk, in the Kondo limit is twofold: First, the position of the 

impurity level is shifted due to the reduction of charge fluctuations, which 

reduces Tk. Secondly, the bare Kondo exchange coupling is enhanced as spin 

fluctuations are enlarged. In total, Tk increases. Both corrections require 

intermediate states beyond the standard Varma-Yafet ansatz. This shows that 

the Hubbard interaction does not just provide quasiparticles, which hybridize 

with the impurity, but also renormalizes the Kondo coupling. 
75.20.Hr, 75.30.Hx, 71.28.+d 71.27.+a 
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I. INTRODUCTION 



Recently, heavy-fermion behavior has been observed in the electron-doped cuprate 
Nd2_a.Ce2,.Cu04 (0.1 ^ x ^ 0.2).i Below 0.3 K a linear specific heat = 7T is observed 
with a large Sommerfeld coefficient 7 ~ 4 J/ (mole Nd ■ K^). In the same temperature regime, 
the spin susceptibility is found to be independent of the temperature and the Sommerfeld- 
Wilson ratio is of order unity. These are characteristic features of heavy-fermion excitations.i 
However, the characteristic low energy scale of the order of 1 K which is associated with 
this behavior cannot be explained by applying the usual theory of the Kondo effect which 
assumes that the conduction carriers behave as free particles.! This is not too surprising 
because undoped Nd2Cu04 is an antiferro magnetic charge-transfer insulator instead of a 
metal,§'i despite of one hole per unit cell. Upon doping the Nd ions are therefore coupling 
to a system of strongly correlated electronsiS rather than to weakly or uncorrelated ones. 

Hamiltonian and Scaling. In order to explain this new type of heavy-fermion behavior, 
it has therefore been proposed to include the correlations among the conduction electrons 
by including an on-site repulsion.lii Thus, the total Hamiltonian 

H = H, + Hf + H,f (1.1) 

goes beyond that of the single-site Anderson impurity model.! is a Hubbard Hamiltonian 
describing the conduction electrons 

He = Ht + Hu 

Ht = ^e(fc)cLcfc^ (1.2) 



k,, 
U 



Hu — X! ■ '^k+S(7'^k(7'^k'-5a''^k'a' 
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creates an electron with spin a and momentum k, is the number of lattice sites. The 
non-interacting dispersion is given by e(A;). : ■■■ : denotes normal ordering with respect 
to the Fermi sea |FS) where all states below the Fermi momentum, kp, are occupied. The 
magnetic impurity is assumed to contain one orbital (e.g., 4/), which is either empty or singly 



occupied. Double occupancies are excluded because of the strong repulsion of electrons in 
that orbital. The energy of the /-orbital is then given by 

Hf = efY.flL, (1-3) 

a 

where = |o")(0| are Hubbard operators forbidding a double occupancy of the impurity 
site and ej < 0. The two subsystems are coupled by a local hybridization 

= E + 4 J.) • (1-4) 

Taking a twofold degeneracy of the /-orbital (cr = 1,2), the model defined in Eqs. 
corresponds! to the situation found in Nd2-xCea;Cu04 since the crystal-field ground 



state of Nd is a doublet .0 In order to perform a systematic expansion we set cr = 1 .. . Nf 
and consider large Nf. This generalization deserves some comment: If the conduction 
electrons are uncorrelated [U = 0), this corresponds to treating an A^/-fold degenerate 
impurity which hybridizes with an s-wave conduction band. This is seen by expanding 
the conduction electron states in partial waves about the impurtity site and, assuming a 
spherically symmetric hybridization, only conduction electrons with the same total angular 
momentum are coupled to the impurity while the others play a passive role and can be 
dropped.0 

Due to the interactions among the conduction electrons this change of basis does not 
simplify the Hamiltonian Nevertheless, we will consider the Hamiltonian in Eqs. (|1 . 1| - 

|1.4|) for cr = 1 . . .Nf, which could be viewed as an SU(A^/) generahzation of the original 
model which has SU(2) symmetry. We thereby create an artificial model, which no longer 
corresponds to the physical situation of an A^j-fold degenerate impurity hybridizing with a 
correlated s-band. The advantage in doing so is, however, that a controlled approximation 
becomes possible for this model, namely an expansion in 

In taking the limit Nf ^ oo, we keep the density of conduction electrons per spin 
constant, so that the kinetic energy increases oc Nf. To have a proper limit Nf oo, the 
hybridization coupling constant V has to be scaled according to V = V/ i/iVj.0 As regards 



the Hubbard interaction, we set [7/2 = U/{2Nf) as was suggested in Ref. |T^. With this 
scahng, the correction to the ground-state energy of the Hubbard model ( |1.2| ) is of order 
Nj, both in second-order perturbation theory in U and when summing the diagrams of the 
random-phase approximation (RPA), which is one order less than the U = energy. 

A straightforward variational ansatz. In case the conduction electrons are uncorrelated 
{U = 0) Varma and Yafet0 proposed the following variational ansatz for the ground-state 
wave function of (|1.1| ) 

l^o)= (l+E«JJv) • (1.5) 

|FS) denotes the filled Fermi sea with empty /-level. Via Hcf, it couples to the states 
/J^Cg^|FS). Each of them describes a singlet formed between the / level and the free electron 
state with momentum q, where q is restricted to occupied states (|g| < k-p). Minimizing 
(^ol-f^— -E'sl^o) with respect to yields an approximate ground-state energy Es- Compared 
to the energy Eyi of the multiplet /J|FS), this collective singlet formation gives rise to a gain 
in kinetic energy. With this energy gain a characteristic temperature scale Tk, the Kondo 
temperature, is associated. In the Kondo limit {\e\ ^ D ^ |ej|, 2D = band width) one 
find^l (in units of the band width) 

TK = nexp(^^) . (1.6) 

Here we assumed a constant density of states p = 1/{2D) of the conduction electrons, fi is 
the chemical potential of the conduction electrons and n = {D + fi)p denotes the filling per 
spin. Subsequently, it has been shown that in an expansion in the inverse degeneracy of the 
magnetic impurity, the ansatz ( p..5| ) yields the ground-state energy to order (1/A'^^)°.E1 

In the case ?7 7^ 0, it is, therefore, tempting to generalize the ansatz ( |1.5| ) by replacing the 
non-interacting ground state |FS) by the (unknown) one of the Hubbard model, \g)^ The 
expectation values with respect to \g) which arise in a variational calculation are given by 
the moments of the spectral function of the Hubbard model, which can be taken from, e.g., 
applying the projection technique.0 If we assume Fermi-liquid behavior for the Hubbard 
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model, the result of this generalized ansatz is obvious. We introduce quasiparticles via 
cj^ = -\/ZcJo- + • • • where Z denotes their renormahzation factor. These quasiparticles 
hybridize with the impurity site rather than bare electrons, the effective hybridization being, 
however, renormalized by \fZ. Therefore, we expect a Kondo temperature 

where pqp is the quasiparticle density of states at the chemical potential. Noting that 
p = ZpQp is the many-particle density of states we see that the correlations enter only via 
p. In particular, Tk is not modified for small U. This is in contrast to Ref. |l^ where it has 
been shown for the Kondo model by a mean-field decoupling that due to polarization effects 
the Kondo temperature increases, even to lowest order in U. In the strongly correlated 
case Eq. ( |1 . 7|) cannot be correct as well since the Kondo exchange coupling should be V'^/U 



rather than V^^/e/.i 

To clarify the quality of the variational approach, we will restrict ourselves to the weakly 
correlated case in this paper and perform a expansion to lowest order in U. The 

theoretical framework, the BriouUin-Wigner perturbation series, is introduced in the next 
section. In particular, we will show that to order l/Nf additional contribution arise from 
the Hubbard interaction in the singlet channel which do not occur in the multiplet channel 
(Sec. PI) and conclude that these contributions modify the Kondo temperature. They are 
estimated in Sec. |^ and the results are discussed in Sec. 0. 



II. BRIOULLIN-WIGNER PERTURBATION THEORY 

The 1/A^j-expansion for the ground-state energy can be derived with the help of Brioullin- 
Wigner perturbation theoryiSlllli We decompose H from Eq. ( |1.1| ) into Hq + Hi and choose 
Hi = Hjj + Hcf as perturbation. In order to obtain the singlet ground-state energy we take 
as unperturbed ground state the filled Fermi sea |FS). The energy of the ground state 
of H (relative to the energy of |FS)) is given bj"^ 
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Es = {FS\H^f:(—^H,] |FS). (2.1) 

Here, Q = 1 — |FS)(FS| and H = Lt + Hf, where the Liouvillean is defined by LfA = 
[Ht,A]^. Equation is equivalent to the zero temperature hmit of the equation for the 
lowest lying pole of the empty /-state propagator that appears in the partition function 



(see, e.g., Ref. |T9[). 

Diagrams. The individual terms of the series (|2.1|) can be visualized by diagrams: In Hi, 
each Hcf changes the occupation of the impurity level from (wiggled line) to 1 (dashed line) 
destroying an conduction electron (solid line), since no double occupancy is allowed (and vice 



versa). This vertex carries a factor V/ ^NfNg. The impurity line changing always between 
occupied and unoccupied /-level constitutes the backbone of a diagram. Hi contains Hu 
as well. The vertex Hu has two incoming and two outgoing conduction electron lines. It 
yields a factor U/ (MfNg) and a ^-function ensuring momentum conservation. Taking the 
expectation value with respect to |FS) we connect the conduction electron lines in all possible 
ways. The resolvent Q/{e — H) yields the energy of the intermediate states and, because 
of the Liouvillean Lt only the energy difference with respect to the filled Fermi sea enters. 
Conduction electron lines pointing to the right correspond to particle-like excitations (with 
a momentum denoted by a capital letter, \Q\ > k-p) while those pointing to the left are 
hole-like (denoted by |g| < kp). Without Hubbard interaction Hu these rules correspond 
to the standard ones0'lli for the self energy of the propagator of the empty state in the 
partition function at zero temperature. 

For an expansion in 1/Nf we note that each closed loop of fermions yields a summation 
over spin and, hence, a factor Nf, whereas each V^- vertex is oc 1/^Nf. To lowest order in 
l/Nf, the application of Hjj does not change the order of a diagram: If we connect two 
conduction electron lines of a closed loop by the 4-point interaction Hu we create two loops 
with the only restriction a ^ a', which is of higher order in 1/Nf. 

Multiplet energy. Similarly to Eq. ( p.l| ), we obtain the energy Eyi for the multiplet 



ground state by taking the expectation values with respect to the multiplet state /j!|FS). 
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The resulting equation corresponds to the lowest pole of the propagator of the occupied 
/-state. The important energy for the low-temperature thermodynamics is given by the 
energy difference of singlet and multiplet ground state, Eu — — -Es, which is related to the 
Kondo temperature.0 Note that the multiplet has one electron more than the singlet in our 
definition. 

Renormalization of the hare propagators. Already to order {1/NfY the bare empty /- 
level (single wiggle line, 1/z) has to be renormalized (double wiggle hne, Go{z)). This 
renormalization arises from a partial summation in Eq. (p.l|) shown in Fig. |l|: 

Go(.) = i + ^/(°)(^)Go(.) = ^3^. (2.2) 
The self-energy J(°) {z) (see also Fig. ^a) evaluates to 

The propagator of the occupied /-state, Gi{z), is not renormalized to this order 

G,{z) = . (2.4) 

z-ef 



III. GROUND-STATE ENERGIES TO ORDER l/Np 

Diagrams for the singlet energy. To order (l/A^/)° only the diagram shown in Fig. 
occurs. It was already evaluated in Eq. 



^^°n^) = ^EG'i(^ + g. (3.1) 

There are no diagrams oc t/ to this order. To order we first find the diagram shown in 
Fig. |6: 

I^'\z) = Y: [Gii^ + e.)]' Go{z + - eg) . (3.2) 

^^f^^s qQ 

As mentioned previously, applying Hu does not change the order of a diagram to lowest 
order in l/Nj. Therefore, I^^^ can be regarded as parent diagram in which we insert vertices 



of the interaction, Hu. Thereby we restrict ourselves to first order in [/, i.e., we apply Hu 
only once in the series We then find the diagrams shown in Fig. ^c. As the diagrams 

are time ordered, J^^"* differs from Iq\ etc. Also, applying Hjj over a doubly wiggled line 
deserves some comment (see, e.g., i'b')'- Such a diagram would not be unambiguous since 
it is not clear whether Hjj acts while the /-level is empty or occupied, when we expand 
the renormalized empty /-propagator as in Fig. ^ For that reason we define that Hu acts 
while the (bare) /-level is empty. The other case yields a different diagram (here, Iq^)- The 
contributions of the diagrams of Fig. ^ are given by 

Ia\z) = - , . , .3 + eq)G'i(z + er + tr' - tR)Go{z + e.^ - tR)Gi{z + e^) 5r'-5,iJ-r 

^^f^^s qrr'R 

Ib\z) = + (^r)Go{z + er- eR)Go{z + - eQ)Gi{z + tq) 6q_Q^r-R 

^^f^^s qrQR 

li^\z) = ^^ E G^iz + er)Go{z + er - eR)Giiz + er + er> - en) 

^^f-^^s qrr'QR 

X Gi{z + eq + er' - eQ)Go{z + - eQ)Gi{z + e^) 5q^Q^r-R 
" E ^o{z + eq + er - eg - eR)Gi{z + + - en) 

^^f^^s qrQR 



X Go{z + er - eR)Gi{z + e^) Sq.q^r-R 



Ie\z) = -^^^ E Gi{z + er')Gi{z + eq + er + er' - eg - eR)Go{z + eq + er - eg - cr) 

^^f-^^s qrr'QR 

X Gi{z + eq + er- eQ)Goiz + - eQ)Gi{z + eg) SQ^q^r-R ■ (3.3) 

According to Eq. ( ^^1]) , the ground-state energy (relative to -Efs) is given by the smallest 
solution of 

E 



Es = I^'\Es) + /(^)(^s) + E ^iiEs) . (3.4) 



i=A 

There is no contribution oc UV^ in this expression for the ground-state energy since we 
introduced the Hubbard interaction in normal ordered form in Eq. ( |1.2D and restricted to 
first order in U. Hence the Hubbard interaction enters only via the hybridization V in the 
ground-state energy. Expanding Eq. (|3.4|) in we obtain 



(0) 



{0)^ 



(3.5) 



l_5/(0)(^(0))/5^g 

Diagrams for the multiplet. We now turn to the ground-state energy of a multiplet state. 
To order (l/Nf)^ it is given by Em = c/ (relative to Epo,)- There is only one diagram 
contributing to order 1/Nf, which is shown in Fig. ^. It is 



(3.6) 



and we find therefore 



EM = ef + J^^\ef) + o{l/Nff . 



(3.7) 



Kondo temperature. We associate the Kondo temperature, Tk, with the difference be- 
tween singlet and multiplet ground-state energyiil (in units of the band width) 



Tk = {Em - fj. - Es)p 



With this definition we find from Eqs. ( p.5|) and ( p.7|) to order (1/A^j)° 



(0) 



(3.8) 



(3.9) 



Assuming a constant density of states we have 



and hence 



D z + e — ef 



pV^ log 



(e/ - P) 



{D + p) + (e/ - p) - z 



Til'' = {ej-p)p-{pVf\o^ 



(0) 



n + Tf 



(3.10) 



where n = {D + p)p denotes the filling per spin. This is solved for small Jk 
by 



(3.11) 

-VViej-p) 
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(0) 



K 



nexp 



(3.12) 



cf. Eq. (|OD. To order we find from Eqs. Q, and Q 



p(l) 

-K 



p 



(3.13) 



l_5/(o)(^(o))/a£;g 

Connection to variational approach. To find the result ( p.l3| ) variationally the following 
states, which occur as intermediate states in the diagrams, have to be included in the trial 
state for the singlet ground state 



^Qc7'^qa'^Qia''^q'a''i ^^^Qcr'^qa'^Q' a''^q' a' fl"'^q" 



(3.14) 
(3.15) 



The variational coefficients are determined up to first order in U and to leading order in 
1/A^j.lli In the free case [U = 0) the first state corresponds to the ansatz of Varma and 
Yafet, cf. Eq. ( |1.5|) , which gives the result correctly to order (l/Nf)^. The next two yield 
the 1/Nf corrections while the last two are of order [l/NjY. 



IV. ESTIMATING THE KONDO TEMPERATURE 

In this section, we estimate the effect of the diagrams oc 1/A^/ on the Kondo temperature. 
We scale the energies by p and study the dependence on T^^ rather than on e'^'' as r^°^ is 
the small quantity. The transformed propagators read 



lyx] 



i-/(°)(Ef -x/p)=iog r^"+^^ 



go{x) = -pV'Go{El^^ -x/p) 



M (4-1) 

x + T^^ - p{ef - p) - p^VH{x) 

and depend implicitly on T^^ and n. The empty-state propagator, go{x) diverges oc /x 

for X — * 0. [This corresponds to the spin fluctuation peak at z = E^^ in Go{z).] This 
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singularity yields contributions oc T^^ (logT^^^^ to and J^^\ which we neglect 

against terms which remain constant as T^^ — > 0. For larger x however, gQ{x) drops slower 
than 1/x resulting in finite contributions. In this intermediate x range (t4°^ < X < n), we 
may safely approximate 

go{x) ~ pJk (4.2) 

for small tI^\ The validity of this replacement for the whole x range in the diagrams has 
been checked numerically. 

Diagrams of order . We begin with the contributions oc f/° in Eq. ( |3.13| ) 



/«(4°)) = -^/ du dxgl{u)g,{u + x) (4.3) 



Nf Jo Jo 

J^'\ej) = --^ dx go{x - Tf ) , (4.4) 

where we again assumed a constant density of states. Inserting (|4.2| ) we find for the multiplet 
energy 

J^'\ef) = -^a-n), (4.5) 

where the corrections are of higher order in l/(e/ — /i)p. For the singlet energy we use the 
same approximation for go to obtain 



Together with the denominator in Eq. ( |3.13|) 



and neglecting the 1, we find that both contributions oc cancel. (Loosely speaking, these 
terms describe the energy gain due to hybridization with unoccupied states which is the 
same for multiplet and singlet state.) 

Diagrams of order U . We continue with the estimation of the diagrams oc U . The 



numerical evaluation of the sums /\ , . . . , J^; is difficult because of the 5-functions, which 



r(i) r(i) 

E 
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ensure momentum conservation in the Hubbard interaction. Since we are interested only 
in the quahtative behavior, we may neglect them. This implies that the interaction U acts 
only at the lattice site with which the impurity hybridizes and corresponds to taking the 
limit of infinite dimensions.^ Then the sums oc f/ in Eq. ( |3.3| ) read 

^2^2 / m(0) 



A (^s ) = — ^y^^ — log ( ^ ^^^(0) J du dx gi{u) go{u + x) i{u + x) 



U r fl—n 1 2 

du / dx gi{u) go{u + x) 



(1) (0) W^p^ ['^-^ 
Ic{El') = — / dudv dx dy gi{u) go{u + x) gi{v) go{v + y) 



(1) (0) 2?7 rl-n 

I}j'{Ey) = — dudv dx dy gi{u) go{u + x) gi{u + V + x) go{u + V + X + y) 



' rn rl—n 

Nf [Jo 

i(u + x) — i(v + y) 
X- — 

u + X — [v + y) 

,^ , du dv 
Nf Jo Jo 

/y(£;g ) = / dudv dx dy gi{u) go{u + x) gi{u + V + x) 

Nf Jo Jo 

X g^[u + + X + ^ ^ . (4.8) 

U + f + X + y 

For the integrals J^^'*, and J^'' it numerically proves sufficient to replace g{x) as in 
Eq. ( [4.2|) in the limit t'^^ — > 0. Then the following leading behavior for t'^^ — > can be 
calculated analytically 

(Es )- j^^ log|^— 

l'^\EP) = ^a-nr (4.9) 

^V/ log 

with C{n) = —2 [nlogn + {1 — n) log(l — n)] > 0. The integral J^^^ simplifies as one inte- 
gration can be performed: 

/W(Ef ) = J_£X!l. rdu dv r^dx g,{u) \ 

^ ^ ^ ' Nf Jo Jo ' i{u + x)-i{0) 

, \ 1 /M + + X + 1 - r2\ ,,^r.S 

x^i M + ^; + x log ■ ■ . 4.10 

V M + f + X / 
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By numerical evaluation, one finds that Iq^ and Ij^^ remain finite as T^""* 0. They are, 
however, small compared to I^^^ as they have an additional prefactor p'^V'^. Keeping only 



diagrams I^^ and I^^ in Eq. (|3.13| ), which corresponds to considering only the states ( p.l4 ) 



in a variational calculation, we obtain in the limit of small Tj^°^ 

4.)._|^U(„).og(f)+^(l-„)= 



(0) 

'^C'H-^(l-n)2) . (4.11) 



The first contribution to T^^ is positive and, therefore, enhances the Kondo temperature. 
Since it depends on f// Jk it is related to spin degrees of freedom. A similar contribution was 
found in Ref. |1^ and it was attributed to the enhancement of spin fluctuations that result 
from the reduction of charge degrees of freedom when turning on U. The second contribution 
in Eq. ([4.11|) depends on pV"^ rather than Jk- It is related with charge degrees of freedom. A 
similar effect has been found in Ref. | and has been interpreted as the increase in energy of 
the virtual state in the spin-exchange process because in the virtual state a conduction site 
is doubly occupied. It decreases the Kondo temperature. However, in the hmit p\ef—p\ ^ 1 
that we considered throughout, the first term dominates: Overall we find an increase of the 
Kondo temperature. 

This interpretation can be put onto more solid grounds by the following observation: If 
we would not scale the Hubbard interaction among the conduction electrons by l/Nj, the 
corrections due to U, I^^\z), would be of the same order as Then the integral I^\z), 
which remains constant as 2; — > 0, would effectively shift the position of the /-level to e*j, 
whereas I^j^\z) ~ logz would renormalize the exchange coupling V'^/e*j [cf. Eqs. ( ^.9[ ) and 



( p.ll| )]. However, without scaling there would be contributions of higher order in U which 
diverge as Nf 00. 



Comparison to previous results. In Ref. |T^, a Kondo model {Nf = 2) with correlated 



conduction electrons has been investigated to lowest order in the interaction strength U by 
a mean-field decoupling of the Kondo-exchange interaction. The following increase of the 
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Kondo temperature has been found 

with a = (3/2)pC/log2. In contrast to our result ( [4.11|) , the increase of the Kondo tem- 
perature seems to depend exponentially on U. Note, however, that in our treatment the 
interaction had to be scaled by 1/Nf. Scaling U in Eq. ( [4 .121) and expanding in yields 
the first term of our result ( [4.11|) (with a factor log 2 for Nf = 2 at half filling instead of 
3/4 log 2). 

The second term of Eq. ( [4.11|) , which describes the effective shift of the /-level, cannot 



be found in Ref. O because there the Kondo model has been investigated, where the charge 



degrees of freedom of the impurity have already been projected out. Therefore, Jk in (|4.12|) 
is an effective coupling constant which depends on f/.i 

Without scaling the Hubbard interaction, all integrals would be of order (1/A^j)°. As 
discussed above, the position of the /-level and the Kondo coupling constant are modified, 
and these corrections occur in the exponent as in Eq. ( [4. 121) . 



V. CONCLUSION 

The aim of this paper was to investigate the influence of correlations among the con- 
duction electrons on the Kondo effect. Lead by the situation prevailing in Nd2-xCea;Cu04, 
we proposed a model with a twofold degenerate impurity {Nf = 2) which hybridizes with 
a correlated s-band and straightforwardly generalized it to arbitrary Nf. As discussed in 
the introduction, this generalization does not correspond to the physical situation of an 
A^j-fold degerate impurity hybridizing with a correlated s-band for Nf > 2, in contrast to 
the uncorrelated case. Although artificial, we saw that this model allows for systematically 
studying the effects of the correlations on those diagrams which are usually considered in 
the uncorrelated case. 

In particular, we assumed that the correlations are weak and calculated their effect on 
the Kondo temperature to lowest order in 1/Nf. We found two competing effects: The first 
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contribution is related to charge fluctuations. Because the energy of the virtual state in 
the spin exchange process increases, the Kondo temperature is reduced. This corresponds 
effectively to a shift of the position of the /-level. A similar effect has been found in Ref. ^. 
The second contribution is related to the enhancement of spin fluctuations of the conduction 
electrons. The Kondo exchange coupling is effectively enhanced and the Kondo temperature 
increases. In the Kondo limit the second contribution dominates the flrst one, so that we 
flnd in total an increase of the Kondo temperature for small U. 

To our opinion more interesting is that corrections to the Kondo temperature occur 
already to lowest order in the Hubbard interaction U. To obtain them in a variational 
approach, trial states are needed which in the uncorrelated case yield corrections to the 
ground-state energy which are of order 1/Nf (and higher). Thus, our result cannot be 
obtained by an ansatz of the Varma-Yafet type (|1.5| ). This shows that the effect of the 
Hubbard correlations is more intricate than just to provide quasiparticles with a modified 
density of states at the Fermi surface which hybridize with the /-orbital as it was described 
by Varma and Yafet for the uncorrelated case, U = 0. 

If we wish to proceed to higher order in U, we note that to order 1/Nf only RPA-type 
diagrams contribute since each U vertex carries a factor 1/Nf, which has to be compensated 
by a spin summation, i.e., a closed loop of conduction electrons (in fact, the ground-state 
energy of the Hubbard model ( |1.2|) to order 1/Nf is given by summing the diagrams of RPA 
type and neglecting the a ^ a' constraint). Only few more intermediate states will occur. 
This is, however, an artifact of our scaling of the Hubbard interaction and one expects that in 
a realistic model (without the restrictive scaling of Hu and finite iVj) intermediate states with 
more and more excited electron-hole pairs contribute, the number of which increases with 
increasing order of U . Therefore, it seems questionable that a systematic treatment 
grasps the correct physics for realistic models of interacting conduction electrons in the limit 
of strong correlations. 
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FIGURES 

FIG. 1. Renormalization of the empty state propagator to order (l/7Vy)° 

FIG. 2. Diagrams for the singlet ground-state energy, a. Order and (l/7Vy)°, I^^\z). b. Or- 
der C/O and il/Nf)\ c. Order and i,l/Nf)\ if (z) {i = A,...,E) 

FIG. 3. Diagram for the multiplet ground-state energy, 
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